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In this work we describe in detail the Chopped RAndom Basis (CRAB) optimal control technique 
recently introduced to optimize t-DMRG simulations [l[. Here we study the efficiency of this control 
technique in optimizing different quantum processes and we show that in the considered cases we 
obtain results equivalent to those obtained via different optimal control methods while using less 
resources. We propose the CRAB optimization as a general and versatile optimal control technique. 
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Realizing artificial, controllable quantum systems has 
represented one of the most promising challenge in 
physics for the last thirty years Q . On one side such sys- 
tems could unveil unexplored features of Nature, when 
employed as universal quantum simulators on the 
other side this technology could be exploited to realize a 
new generation of extremely powerful devices, like quan- 
tum computers Q. Along with the impressive progress 
marked recently in the construction of tunable quantum 
systems [|J @, there is a renewed and increasing inter- 
est in quantum optimal control (OC) theory, the study 
of the optimization techniques aimed at improving the 
outcome of a quantum process 0- Indeed OC can prove 
to be crucial under several respects for the development 
of quantum devices: first, it can be generally employed 
to speed up a quantum process to make it less prone to 
decoherence or noise effects induced by the unavoidable 
interaction with the external environment. Second, con- 
sidering a realistic experimental setup in which just few 
parameters are tunable or, in the most difficult situa- 
tions, only partially tunable, OC can provide an answer 
about the optimal use of the available resources. 

Traditionally OC has been exploited in atomic and 
molecular physics 0-0- More recently, with the advent 
of quantum information, the requirement of accurate con- 
trol of quantum systems has become unavoidable to build 
quantum information processors [Tol - fT6| . However, the 
above mentioned methods often result in optimal driv- 
ing fields that require a level of tunability incompatible 
with current experimental capabilities and in general, the 
calculation of the optimal fields requires an exact descrip- 
tion of the system (either analytical or numerical). The 
field of application of these methods is severely limited 
also by the need to have access to huge amount of infor- 
mation about the system, e.g. computing gradients of the 
control fields, expectation values of observables as a func- 
tion of time. Moreover, standard OC algorithms define a 
set of Eulcr-Lagrangc equations that have to be solved to 
find the optimal control pulse Q , where the equation for 
the correction to the driving field is highly dependent on 
the constraints imposed on the system and on the figures 
of merit considered. This implies that considering dif- 
ferent figures of merit and/or constraints on the system 
needs a redefinition of the corresponding Eulcr-Lagrangc 
equations, hindering a straightforward adaptation of the 



optimization procedure to different situations. 

In this work we discuss in detail the Chopped RAn- 
dom Basis (CRAB) technique, an optimization method 
directed to overcome these difficulties and already in- 
troduced in [l|. The CRAB optimization is based on 
the definition of a truncated randomized basis of func- 
tions for the control fields that recast the problem from 
a functional minimization to a multi-variable function 
minimization that can be performed, for example, via a 
direct-search method. As shown in the following, the 
CRAB optimization flexibility allows to construct OC 
pulses just exploiting the available resources. Indeed, 
different figures of merit and constraints can be easily 
considered without any complications. Another appeal- 
ing characteristic of CRAB is its compatibility with t- 
DMRG techniques: this feature indeed significantly en- 
larges the class of controllable systems [l[ , from few-body 
or exactly solvable to general many-body quantum sys- 
tems with "moderate" deg ree of entanglement generated 
during the dynamics [l7|- This is, to the best of our 
knowledge, the unique OC algorithm that can be applied 
in such vast setting. Finally, it can be straightforward 
applied also in a closed- loop optimization experiment, 
where the simulation of the system under study is re- 
placed with the experiments itself. 

Here we analyze the CRAB optimization as a possi- 
ble general OC algorithm to be used also in a standard 
context (solvable and/or few body systems) as a valid 
alternative tool with respect to standard OC methods 
to find optimal control fields. Indeed, recently optimiza- 
tion methods based on the expansion over a particular 
function basis have shown to be effective (l8N2l| . In par- 
ticular, a similar approach has been proved to be math- 
ematically convergent and consistent [22|, [23| . On top of 
that, some theoretical analysis over control landscapes 
suggests that, at least in the absence of constraints, the 
figure of merit landscape might be smooth enou gh t o al- 
low for simple optimization procedures to work |24l. |25| . 
Here, we show that indeed a convenient choice of the 
function basis driven by physical or geometrical argu- 
ments is enough to obtain optimal driving fields. How- 
ever, in the cases where no physical intuition drives the 
choice of the function basis, the CRAB algorithm allows 
to find the optimal driving fields where a simple ansatz 
would fail. Moreover, a comparison between the results 
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of CRAB with and without a physically driven choice of 
the basis, as well as previous results obtained using dif- 
ferent optimal control algorithms (Krotov's algorithm), 
show comparable performances (26[. 

The structure of the paper is the following: in Sec.Uthe 
CRAB optimization is described; in Sec. [II] it is applied 
to a paradigmatic quantum control problem, the state 
transformation of two coupled qubits, to show its poten- 
tial. Then we compare the results obtained via CRAB 
optimization in more complex cases already present in 
literature [26|, H?}: in Sec. IIIII the method is employed 
to control the quantum phase transition evolution of the 
Lipkin-Mcshkov-Glick (LMG) model; and in Sec. IIVI we 
optimize the transfer of a state along a spin chain. Fi- 
nally, in Sec. |V] the optimization is exploited to maxi- 
mize the final entanglement entropy of the final state in 
the LMG model; and in Sec. I VII a comparison between 
adiabatic and optimized processes is proposed. 



I. CRAB OPTIMIZATION 

The optimization problem we are dealing with is de- 
fined as follows: given a Hamiltonian H acting on a 
Hilbcrt space H = C , depending on a set of time- 
dependent driving fields T(t), we search for the opti- 
mal transformation to drive, in time T, an initial state 
V'o) £ H into a different one (target state) \i/>c) £ H 
with some desired properties expressed by a cost func- 
tion /(|-0g)) we want to minimize [55j . In addition, con- 
straints might be present on the driving fields, e.g. to 
match experimental conditions: They can be expressed 
usually as a function of the driving fields Ci(T(t)). Typ- 
ical scenarios and corresponding cost functions and con- 
straints are: 

1. The goal is the preparation of a well-defined quan- 
tum state IV'g) with high accuracy for which a con- 
venient cost function is the infidelity, 

hmT)))^i{T) = i-mT)\^ G )\\ (i) 

2. The target state is the unknown ground state of a 
Hamiltonian H p . The cost function is then given 
by the final system energy, 

/ 2 (|V(T))) = Ef (T) = (ip(T)\Hp\ip(T)) . (2) 

3. The target is some property or condition that many 
states can satisfy, like for example, in the produc- 
tion of highly entangled states. In this case the cost 
function is simply defined as 

/ 3 (|V(T))) = -S(MT))), (3) 

where S(\^}) is a convenient measure of the entan- 
glement of the state 



4. A constraint is present on the power of the driving 
fields, that is, the solution should minimize also the 
fluences 




5. A limited bandwidth is allowed for the driving 
fields: below we show how this is already embedded 
in the algorithm and is not necessary to consider it 
as an additional explicit constraint. 

6. The initial state or the driving fields are known 
within a given uncertainty e. In this case, the cost 
function can be defined as an average other all pos- 
sible outcomes compatible with that uncertainty, as 
for example: 

h = J f(MT,e)))<k. (5) 

All of the aforementioned optimization problems are 
then recast in the problem of solving the Schbdingcr 
equation (from now on we assume ft = I) 

i^ t m)) = H[r(t)}\m), (6) 

with boundary condition \ipi) — |V>(0)), while minimizing 
the cost function 

F = af + J2&C l (r(t)) 7 (7) 

i 

where the coefficients a and Pi allow for a proper weight- 
ing of the different contributions (the /3s play the role of 
Lagrange multipliers) and / is the chosen cost function. 

To perform such an optimization, the CRAB algorithm 
starts from an initial pulse guess r°(£) and then looks for 
the best correction of the form 

r CRAB (t)=r (t) . 5j(t) , (g) 

The functions gj(t) are expanded in a simple form in 
some function basis characterized by some parameters 
VLj (Fourier space, Lagrange polynomials, etc.): gj — 
Y, k c|<$(fif). The two key ingredients of the CRAB 
optimization are that the function space is truncated to 
some finite number of components N c (k = l,...,N c ) 
and that the corresponding basis functions are "random- 
ized" to enhance the algorithm convergence, i.e. — > 
gj(flj(l+rj)) where is a random number. Indeed, this 
last choice breaks the orthonormalization of the functions 
gj, however as we show in the following, it allows for an 
improved convergence of the algorithm as it enlarge the 
subspace of functions explored by the algorithm while 
keeping constant the number of optimization parameters. 

The optimization problem is then reformulated as the 
cxtremization of the multivariable cost function J-(T, 5j), 
which can be numerically approached with a suitable 
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FIG. 1: Infidelity /i of the final state as a function of num- 
ber of calls to the optimization algorithm Nj for two capaci- 
tively Josephson charge qubits with principal harmonics (dark 
grey [blue] line) and randomized frequencies (light grey [green] 
line), for the goal state \ipc) an d N c = 2 for thirty different 
random instances. 
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FIG. 2: Optimized infidelity / as a function of the number of 
optimization parameters N c with principal harmonics (dark 
grey [blue], full symbols) and randomized frequencies (light 
grey [green], empty symbols) for different target states \4>g) 
(circle), \ipa) (squares), \if>a) (diamonds). 



method, e.g., steepest descent, conjugate gradient or di- 
rect search methods [28|. Hereafter we use the last op- 
tion, which is the simplest one and easily compatible with 
any technique employed to solve the dynamics induced by 
if[r(t)] (either exact solution of the Eq. © or approx- 
imate solution with time dependent DMRG [ItJ )• This 
choice also gives another advantage with respect to other 
OC methods where gradients and functional derivatives 
have to be computed, increasing the complexity of the 
optimization procedure. 

As an example, in the following problems, we focus on 
the case of a single control parameter T(t) and we choose 
to work in the Fourier basis. The optimal pulse can then 
be written as 

[J2n=i A n sm(uj n t) + B n cos(w„i) ) 



where X(t) is a time dependent function enforcing the 
boundary conditions (i.e. X(t) — > oo for t — > and for 
t -¥ T). The function r CRAB (i) is fixed by selecting 
the optimization parameters A, B and uj, with N c the 
dimension of each vector. In conclusion, given a fixed 
total evolution time T, the cost function is clearly just a 
function of the control parameters, 

T = T crab {A,B,lo). (10) 

The optimization problem is reduced to the minimiza- 
tion of 7 rCRAB (A, B : us) as a function of 3 x N c vari- 
ables. As mentioned before, however, the space of the 
variables can be reduced even more: although in princi- 
ple the frequencies uj can be considered free variables it is 
often convenient to keep them fixed and to perform the 
minimization just with respect to A and B. Indeed as 



shown in our analysis this is sufficient to obtain good re- 
sults. In this approach we need then a criterion to select 
the w's. When we have no available information about 
the typical energy scales of the system under consider- 
ation, the frequencies are picked randomly around prin- 
cipal harmonics: uik = 27rfc(l + r^/T, with random 
numbers with flat distribution in the interval [—0.5,0.5] 
and k = 1, N c . Viceversa when the physical details of 
the model are known, clearly one can exploit this infor- 
mation to select the relevant frequencies, as shown in the 
following sections. 

II. TWO-QUBITS OPTIMIZATION 

In this section we apply the CRAB optimization to 
a paradigmatic problem in quantum information theory 
and control: we search for the optimal way to perform a 
state transformation of a two-qubit system, in particular 
we consider two capacitively coupled Josephson charge 
qubits, even though the following analysis can be easily 
adapted to different qubit implementations. The Hamil- 
tonian of the i-th qubit is defined as [29l . [30| 

H t = E c a\ + Ejal 

where the as are Pauli matrices, Ec is the charging en- 
ergy and Ej is the Josephson energy and i = 1,2. For 
capacitivc coupled qubits, the interaction Hamiltonian 
reads 

Hi = E cc o- 1 z <j1, 

where E cc is the charging energy associated to the 
Coulomb interaction between the qubits. Hereafter we 
set Ej/Ec = — 1, while the coupling will be the driving 
field E cc {t)/Ec = T(i) we use to optimize the transfor- 
mation. We will consider as initial state the state with 
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FIG. 3: (Color online) Infidelity as a function of the size in 
the LMG model. Squares represents the data before the opti- 
mization, circles the data after the optimization with CRAB. 



no excess Cooper pairs | Vo) = 1 00), and our goal states 
will be three different state with different properties: the 
reversed separable state \ipo) = [11), the homogeneous 
superposition state = ^Ylijlhj)' anc ^ the maxi- 

mally entangled Bell state \^%) = ^(|00) + |11)). Note 
that due to the fact that only the coupling is controlled, 
all three states are not trivial to achieve. We set the to- 
tal time of the transformation to the somehow arbitrary 
time scale T = tt/Ej and we perform a CRAB optimiza- 
tion using the truncated expansion of the function g(t) 
given in Eq. Q, with a constant initial guess for the driv- 
ing field r (t) = r(0) = 1. We considered an additional 
constraint on the fluence of the control field, thus the 
resulting cost function is defined as 



^ = /i + 0.1Ci(T(t)), 



(11) 



where f\ and C\ are given by equations ([T]) and (U]) re- 
spectively. Here we are interested in studying the effect 
of the randomness introduced in the frequencies of the 
expansion (jH]), thus we optimize both in the case of ran- 
dom ffc and with ru = 0. To perform a fair compari- 
son, we ran the optimization in both cases with the same 
maximum number of calls Nf ~ 30.000 to the function 
J 7 , which fixes the simulation complexity. Indeed, in the 
first case we repeated the optimization for thirty different 
Tk random configurations (with a single Ak , Bk random 
starting point), while in the second case the optimization 
was repeated over thirty initial random Ak, Bk configura- 
tions. A typical result is shown in Fig. [TJ for N c = 2 and 
\i/jq): it clearly shows that for the case of randomized 
u>k the optimization is highly improved (notice the loga- 
rithmic scale). A more systematic comparison is shown 
in Fig. [2] where the best results are plotted against the 
number of optimization parameters N c for the three tar- 
get states \iPq): in all cases the randomization of the 
frequencies improves the convergences to higher fidelities 
up to the simulation error. In particular, in one case, the 
final result without randomization is very far from being 




FIG. 4: (Color online) Infidelity as a function of the number of 
control parameters for different sizes in the LMG model. The 
total evolution time is T = 2Tqsl = 2tt/A. Inset: infidelity 
as a function of the number of parameters for a single size 
N = 32: comparison between data optimized using as cost 
function the infidelity (empty circles) and the final energy (full 
circles). Green squares represent the results with randomized 
frequencies. 



satisfactory as the final fidelity is of the order of ten per- 
cent, resulting in a very poor state transformation. On 
the contrary, using the randomized frequencies we were 
able to find optimal pulses to obtain fidelities below one 
percent - values that are comparable, in most cases, with 
experimental errors. 



III. LIPKIN-MESHKOV-GLICK MODEL 

The Lipkin-Meshkov-Glick (LMG) model is the 
paradigm of a system with long range interaction (in- 
finite in the thermodynamical limit). The Hamiltonian 
is written as [3lL S3 : 



H 



-E 

TV ^ 

i<j 



N 

r(t)5>? 



(12) 



where J is the uniform spin-spin interaction (we set ,7=1 
in the following) , N is the number of spins in the system, 
r is the transverse field and erf are the Pauli matrices. 
By introducing the total spin operator S a = ^2i<yf /2, 
Eq. (|12p can be rewritten, apart from an additive con- 
stant, as H = -jf[Sl + yS%] - TS Z . The Hamilto- 
nian hence commutes with S 2 and does not couple states 
having a different parity in the number of spins point- 
ing in the magnetic field direction: [H, S 2 ) = and 
[H, Y\ i of] =0. In the isotropic case 7 = 1, also the 

z-component of S is conserved, [H,S Z ] = 0. In the ther- 
modynamical limit the LMG model undergoes a second 
order quantum phase transition at r c = 1 from a para- 
magnet (r > 1) to a fcrromagnct (r < 1). The phase 
transition is characterized by mean-field critical expo- 
nents [32j|. The phase transitions dramatically affects 
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FIG. 5: (Color online) Infidelity as a function of the size in 
the transfer state problem. Squares represents the data before 
the optimization, circles the data after the optimization with 
CRAB. 



the dynamical behavior of quantum systems: As dis- 
cussed in more detail in Sec. I VII the gap closure at the 
critical point promotes dynamical excitations, preventing 
adiabatic evolutions whenever the adiabaticity condition 
T ;§> A -1 is not fulfilled, where T is the total evolution 
time and A the minimum spectral gap [33TJ42T ] . Following 
Ref. [26j], we employ the CRAB optimization to drasti- 
cally reduce the residual density of defects present in the 
system in a strongly non adiabatic dynamics, drastically 
reducing the time needed to connect the ground state 
in one phase with the ground state of the other phase 
with respect to adiabatic non-optimized strategics. We 
chose as initial state the ground state (gs) of H[T(t)} at 
r, 3> 1, i.e. the state in which all the spins are polar- 
ized along the positive z-axis (paramagnetic phase). As 
target state we chose the gs of H\T = 0] (ferromagnetic 
phase). We focused our attention on the case 7 = 0, 
representative of the class 7 < 1 (for 7 = 1 the dy- 
namics is trivial due to the symmetry of H) [43| . For 
this model indeed a lot of physical information is avail- 
able: the gap between the ground state and the first ex- 
cited state closes polynomially with the size at the crit- 
ical point [32], A ~ N~ x / 3 . Furthermore it has been 
recently demonstrated that the minimum time required 
to obtain a perfect conversion between the initial and the 
final state here considered, the so called quantum speed 
limit, is given by Xqsl = 7r/A [26|, [44j. In order to test 
the performance of CRAB, we fixed the total evolution 
time above this threshold, at T = 2Tqsl, in a regime in 
which in principle it is possible to produce an arbitrarily 
small infidelity with optimized evolutions. 

The results of our simulations for the LMG model are 
summarized in Fig. [3] and Fig. 2] the data shown in the 
two pictures (with the only exception of the inset of Fig. @] 
as explained in the following) have been produced as- 
suming Eq. © as control field and the infidelity as cost 
function to minimize. In Fig. [3] we plotted the infidelity 
as a function of the size N, before the optimization for 
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FIG. 6: (Color online) Infidelity as a function of the number 
of control parameters for different sizes in the in the state 
transfer model. The total evolution time T = 2Tqsl- Inset: 
infidelity as a function of the number of parameters divided 
by the size. 



a linear driving field T°(t) oc t/T (squares), and after 
the optimization with CRAB (circles): for each size we 
have been able to produce an infidelity below 10~ 6 start- 
ing from an infidelity of order 0(1). In particular the 
data have been produced by minimizing Eq. (|10[) with 
respect to A and B, while keeping u fixed, for a total of 
2 x N c = 16 parameters. In this case the frequencies Co 
have been chosen by exploiting the physical information 
available. We chose the frequencies equal to the mini- 
mum spectral gap u>i = 2ir/T = 27t/2Tqsl = A and we 
considered the main harmonics uj^ = kuii for k up to 
N c . In Fig. @] we plot the infidelity as a function of num- 
ber of parameters employed to build the optimal field of 
Eq. ([9]) - adding a frequency ujk corresponds to add two 
parameters, Ak and B^. First it can be noticed that 5 
harmonics are sufficient to reach the best optimization 
result, I ~ 10~ 6 ; however with only 3 harmonics the 
infidelity is already of order 10 -4 , of the order of the 
required threshold for fault-tolerant quantum computa- 
tion. Considering the implementation of an optimal pulse 
in an NMR or quantum optics experiment, the gain with 
respect to other OC methods providing a totally arbi- 
trary r opt (i) is evident. The second interesting feature 
is that the behavior of the infidelity in Fig. @]is approx- 
imately independent of the size (for the smallest system 
considered, N = 10, finite size effects are more evident): 
this confirms the intuition that the most relevant energy 
scale for the LMG model is given by the minimum spec- 
tral gap. 

Finally, in order to verify the independence of the op- 
timization from the knowledge of the target state, we 
repeated the simulations assuming as a cost function the 
final energy E f (T) of Eq. ©. In the inset of Fig. H 
we compare the infidelity of the data optimized using as 
cost function the infidelity itself (empty circles) and the 
final energy (full circles) , for a specific size of the system 
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FIG. 7: (Color online) Final entanglement entropy as a func- 
tion of the total evolution time T for iV = 10, 32, 100 of one 
of two equal bipartitions of the system. The time is measured 
in units of J~ . 

N = 32 and for different number of control parameters: 
as shown in the picture the agreement is very good. We 
also repeated the optimization using randomized frequen- 
cies, obtaining the same results as before. Thus, also in 
the case where the chosen frequencies are optimal, in- 
troducing randomness does not prevent the optimization 
to work. On the contrary, if one has no access to any 
information on the system, the randomization does not 
prevent to reach the same optimal result. 

IV. STATE TRANSFER ALONG A SPIN CHAIN 

In this section we study the optimization of a model 
representing a possible implementation of a quantum bus. 
The model consists in a chain of spins coupled via uni- 
form nearest-neighbor (n.n.) interaction; by acting with 
an external, parabolic magnetic field it is possible to 
transfer a quantum state along the chai n fToL l27l |45| . 
In particular we follow the lines of Ref. (45l. |46| . The 
Hamiltonian of the system is 

T N-l N 

H W = ~ g E °» • ^ + E B n(t)<, (13) 
n— 1 n—1 

where N is the number of spins in the chain, <r„ repre- 
sents the Pauli n— th-spin operator, J is the uniform n.n. 
interaction (we set J = 1 in our simulations), and B n (t) 
is the tunable magnetic field along the z-direction. In 
particular we considered a parabolic magnetic field tun- 
able in position and strength [45j . 

B n (t) = C(t)(x n -d(t)) 2 , (14) 

where d(t) is the position of the potential minimum along 
the chain, x n is the position of the n— th spin, and C(t) 
is the instantaneous curvature of the field. Far from the 
minimum, the spins are forced by the magnetic field to be 



FIG. 8: (Color online) Entanglement entropy saturation 
value as a function of the size (red circles) and the function 
,41og 2 (iV/2 + 1) (dashed line). A fit gives A = 0.947. 



aligned along the z-axis irrespective of their mutual inter- 
action; instead close to the minimum, the n.n. coupling 
prevails and can be exploited to transfer the information 
(i.e. the state) from one site to the next. The Hamil- 
tonian commutes with the total magnetic field along the 

z-direction, [H{t),Y^n=i a n\ ~ ®> so tnat tue dynamics 
occurs in a subspace whose dimension grows just linearly 
with the size N of the system. We chose to work in 
the subspace (J2n=i a n) = li m particular we aimed at 
transferring a spin-up state from one end of the chain 
to the opposite end, or in other words to transform the 
state = |10...0) into the state \ip G ) = |0...01), with 
(1) corresponding to the nth spin pointing in the down 
(up) direction along the z-axis. We employed CRAB to 
optimize the two control parameters, Ti(t) = d(t) and 
T 2 (£) = C(t); as in the previous section, we set the total 
evolution time above the quantum speed limit threshold 
at the value T = 2Tqsl, where for the latter we used 
the estimate made in Refs [13, H||. The optimization 
has been performed by keeping C3i,uj2 fixed (in particu- 
lar LUik — L02k = 2kn/T for k = 1, N c ) and minimizing 
the infidelity with respect to Ai,B ll A 2 ,B2 1 where the 
index 1 and 2 refer to d(t) and C(t) respectively. 

The results of our simulations for the state transfer 
along the chain are summarized in Fig. [S] and Fig. [S] In 
Fig. O we show the infidelity as a function of the size 
before the optimization (squares), for a constant C(t) 
and d(t) = t/T, and after the optimization with CRAB 
(circles): for each size considered we were able to reach 
an infidelity below the value 10 -4 starting from an ini- 
tial infidelity of order 1. In Fig. [6] we plot the infidelity 
as a function of the number of parameters employed in 
the minimization procedure; in this case, unlike for the 
LMG model in Fig 01 the data show a strong dependence 
on the size. We interpreted this behavior as a conse- 
quence of the structure of the problem. Considering the 
particular transfer mechanism, in which the information 
moves step by step from one site to the next one, we ex- 
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pect the optimal pulse to be able to modulate the mag- 
netic field around each spin; this occurs only when the 
spectrum of the pulse involves frequencies of the order 
of the inverse of the time spent on a generic site n, i.e. 
lj ~ 2n/(T/N) = Nun. As a test, in the inset of FigH 
we plotted the infidelity as a function of the number of 
parameters divided by the size; the good agreement of 
the rescaled data confirms our expectation. 



V. ENTANGLEMENT ENTROPY 
MAXIMIZATION 

Among its various applications, OC can be exploited 
for entanglement production pTl . |48| . Here we employ 
the CRAB technique in the LMG model to maximize the 
von Neumann entropy Sl,n = — Tt(pl,jv log2 Pl.n) as- 
sociated to the reduced density matrix pl,n of a block 
of L spins out of the total number N at a given fi- 
nal time T, which gives a measure of the entanglement 
present between two bipartitions of a quantum systems. 
As seen in Sec. IIIII due to the symmetry of the Hamilto- 
nian [H, S 2 } = 0, the dynamics is restricted to subspaces 
with fixed total angular momentum; in particular assum- 
ing as initial state the ground state of the system, we 
have S = N/2. The Dicke states \S = N/2,S Z ) with 
S z = —N/2,..., N/2 provide a convenient basis span- 
ning the subspace accessible through the dynamics. In- 
deed the entanglement entropy Sl,n can be easily eval- 
uated noticing that, since the maximum value of the to- 
tal spin can be achieved only with maximum value of 
the spin in each bipartition, the following decomposition 
holds [H,|i| : 

L 

\N/2,n) = Y,ptg\L/2,l-L/2) (15) 
z=o 

® \(N-L)/2,n-l-(N-L)/2), 




t 



FIG. 9: (Color online) Instantaneous excitation probabilities 
Pi of the i— th excited level (Pi continuous, P2 dotted, P3 dot- 
dashed, P4 dot-dot-dashed, P25 dash-dash-dotted line) and 
total excitation probability P tot = YliLi P (dashed purple 
line) in the LMG model with N = 50 for an evolution induced 
with a driving field linear in time, T(t) oc —t/T, T/Tqsl = 2, 
N = 50. The thick red dashed line signals the crossing of the 
critical point. The time is measured in units of J -1 . 



a behavior closely resembles the features observed in one- 
dimensional systems after a sudden quench [50j | , although 
here we are dealing with a fully connected model [5l|, ■ 
In Fig. [8] we plotted the saturation value reached with 
the optimization as a function of the size N; comparing 
our data with the maximum possible value of the von 
Neumann entropy for a subsystem of L = N/2 spins 
(described by a Hilbcrt space of dimension N/2 + 1) 
= log 2 (A/2 + 1), we obtain almost the maximal 
possible amount of entanglement, S op t/ S max ~ 0.95. 



VI. LINEAR VS OPTIMAL DRIVING 



where n and I correspond respectively to the number of 
up spins in the whole system and in the block of size L, 
andpj,„ = Ll(N - L)\n\{N -n)\/[l\(L-l)\(n-l)\(N L - 
n + iy.Nl}. Expressing the evolved state \ip(T)) in the 
Dicke state basis and using the previous decomposition, 
it is immediate to evaluate Sl,n(T). 
In our simulations we considered a system equally bi- 
partite, i.e. L = N/2, and we took as starting state 
the ground state of the LMG Hamiltonian at T ^ 1, in 
which all the spins are polarized along the positive z di- 
rection, so that the state factorizes and the entanglement 
entropy vanishes, see Fig. [7] Then we performed the op- 
timization with CRAB, modulating the field according 
to Eq. © and using as a cost function Eq. ([3]). The 
behavior of entanglement entropy after the optimization 
S pt(T) for different values of the total evolution time 
T is shown in Fig. [7J after a short transient of linear 
growth, S op t (T) reaches a saturation value growing with 
the size, as expected. It is interesting to notice that such 



In this section we analyze in more detail the features 
characterizing the optimal dynamics induced by CRAB. 
In order to better understand the matter, we draw a 
comparison with a simpler non-optimized dynamics, in 
which the driving field is linearly dependent on time; in 
particular we focus the attention on the LMG model. 
An important point in the study of the dynamics of a 
quantum system is usually represented by the adiabatic 
theorem |33|. The latter establishes that a system ini- 
tially prepared in its ground state can be driven by a 
time dependent Hamiltonian adiabatically (i.e. without 
introducing excitations), if the time scale of the evolu- 
tion is much larger than the minimum spectral gap, i.e. 
T 3> A -1 . In critical systems the spectral gap closes 
at the phase transition, so that the system gets excited 
from the instantaneous gs while crossing the critical point 
for any finite-time evolution [34| . For finite-size systems, 
the critical gap is not completely closed, but it presents 
a pronounced minimum where the excitation appears, 
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as shown in Fig [9] an estimate of the excitations in- 
duced by a linear driving can be obtained by Kibble- 
Zurek theory [3fl [13, EH 0, EH . In the picture we mon- 
itored the instantaneous total excitation probability Ptot 
(dashed line), and the populations of lowest levels (differ- 
ent style [color] lines) during the dynamics. The evolu- 
tion starts at large negative times (left) and ends at the 
time t = (right) ; the critical point is crossed around the 
time t = 11 when T(t) ~ 1, see section Unl Far from the 
critical point the system evolves adiabatically as demon- 
strated by the low total instantaneous excitation proba- 
bility; notice that before reaching the critical point the 
total excitation probability coincides with the small exci- 
tation of only the first level (red continuous line) . In a re- 
stricted region around the critical point (—15 < t < —10) 
the total excitation probability jumps to values of order 1 
and does not change significantly any more. Notice that 
in the final part of the evolution more levels get popu- 
lated, as shown by the difference between the instanta- 
neous infidelity and the excitation probability of the first 
level. At the final time t = T the excitation probability 
is equal to the infidelity of the process, i.e. Ptot{T) = f\. 
We then optimize the final infidelity, and the correspon- 
dent plot for the optimal evolution is reported in Fig [TO] 
The scenario in this case is completely different: the sys- 
tem is excited at the very beginning of the dynamics and 
remains excited for the most part of the evolution until 
close to the end, when the infidelity drops abruptly to 
zero. It is interesting to notice that just a few levels are 
excited, as demonstrated by the small difference between 
the total excitation probability (dashed line) and the ex- 
citation probability of the first level (red continuous line) . 
This result is in agreement with previous findings where 
the authors showed that this kind of dynamics can be 
approximated by a two-level system dynamics J2f|. The 
abrupt jump in the probabilities around the time 20 is 
due to an abrupt (double) change of sign in Fcrab (i), re- 
versing suddenly the order of the levels and transforming 
the gs in the most excited state (dash-dash-dotted [cyan] 
line) , thus this signature is not due to a collective involve- 
ment of all the levels but simply to a reshuffling of their 
order. Indeed as shown in the picture for —18 < t < 0, 
with the subsequent change of sign the previous order 
is reestablished. We can then summarize the main fea- 
tures of the optimal evolution induced by CRAB in three 
points: it is strongly non adiabatic; it involves just a re- 
stricted number of levels although not necessarily close 
to the nominal instantaneous ground state and it is such 
that at the very end all populations constructively inter- 
fere to obtain the desired goal state. 



VII. CONCLUSIONS 

In this paper we studied in detail the performance of 
quantum optimal control through the CRAB optimiza- 
tion [l[. In particular we focused the attention on three 
different systems and different figures of merit, in order 




FIG. 10: (Color online) Optimal instantaneous excitation 
probabilities Pi in the LMG model for T/Tq S l = 2,N = 50. 
Codes are the same as for Fig. [9] The time is measured in 
units of J -1 . 



to outline the versatility of the method. We first studied 
the optimization of state transformations of two qubits 
via a controlled coupling. We have shown that the CRAB 
optimization is very effective already using only a few op- 
timization parameters and the fundamental role that the 
randomization of the function basis plays in increasing 
the process convergence. We then analyzed two many- 
body quantum systems: the first one, the LMG model, 
is the prototype of many-body system with long range 
interaction undergoing a quantum phase transition. The 
success of CRAB in this context confirms the possibility 
of controlling complex systems typically studied in con- 
densed matter, with relatively small resources: due to 
CRAB unique features, only few parameters (3 frequen- 
cies) are indeed sufficient to obtain excellent results. The 
second many-body quantum system studied, the transfer 
of information along a spin chain, is a typical problem 
studied in quantum information theory: the high accu- 
racy achievable through CRAB optimization makes it a 
valuable tool for this kind of applications. Moreover, due 
to the simple structure of the optimal pulses, they may 
be used to extract information on the typical timescales 
involved on the system dynamics, as we did for the infor- 
mation transfer in spin chains. We stress also that the 
exponential dependence of the figures of merit as a func- 
tion of the number of parameters found in all cases (see 
Figs. [21 mill) suggests that in general already a moderate 
number of optimization parameters will be sufficient to 
get huge improvements in the desired processes. 

Finally, we have shown that with a simple change of 
the cost function, the CRAB optimization can be used 
to optimize the search of the unknown ground state of a 
Hamiltonian or to generate quantum states satisfying de- 
sired properties, i.e. high entangled states. Monitoring 
the instantaneous excitation probabilities generated by 
the optimized process, we have demonstrated the highly 
non adiabatic character of the dynamics and the fact 



9 



that, despite the complexity of the system under study, 
just a restricted number of excited levels are really pop- 
ulated during the evolution. The latter fact justifies the 
compatibility of CRAB with DMRG-like techniques. We 
mention that the CRAB optimization has been applied 
also to open quantum systems obtaining interesting re- 
sults and thus increasing its possible applications |54j . 

In conclusion, the main features of the CRAB opti- 
mization -versatility (different constraints, compatibility 
with approximate simulation methods and experiments) , 
fast convergence (the final error scales exponentially with 
the number of optimization parameters while the num- 



ber of algorithm iterations linearly) and simplicity (small 
modification to existing numerical codes for quantum sys- 
tem simulations)- demonstrate that the CRAB optimiza- 
tion is not only an unique solution for many body quan- 
tum systems optimal control but it is a valid alternative 
also in many different settings where other optimal con- 
trol tools exist @. 

We acknowledge discussions with R. Fazio, V. Giovan- 
netti, M. Murphy, and G. Santoro, and support from the 
EU projects AQUTE, PICC, the SFB/TRR21 and the 
BWgrid for computational resources. 



P. Doria, T. Calarco, and S. Montangero Phys. Rev. Lett. [25 
106, 190501 (2011). 

C. Brif, R. Chakrabarti, and H. Rabitz, New J. Phys. 12, [26 
075008 (2010). 

R. Feynman, Int. J. Theor. Phys. 21, 467 (1982). [27 
S. Lloyd, Nature 406, 1047 (2000). 

I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). [28 

G. Roati, C. D'Errico, L. Fallani, M. Fattori, C. Fort, 
M. Zaccanti, G. Modugno, M. Modugno, and M. Ingus- 

cio, Nature 453, 895 (2008). [29 
A. P. Peirce, M. A. Dahleh, and H. Rabitz, Phys. Rev. 
A 37, 4950 (1988). [30 
T. Calarco, U. Dorner, P. Julienne, C. J. Williams, and 
P. Zoller, Phys. Rev. A 70, 012306 (2004). 
N. Khaneja, T. Reiss, C. Kehlet, T. Schulte-Herbruggen, 
and S. G. Glaser, J. Magn. Res. 172, 296 (2005). [31 
S. Bose, Phys. Rev. Lett. 91, 207901 (2003). 
S. Montangero, T. Calarco, and R. Fazio, Phys. Rev. [32 
Lett. 99, 170501 (2007). [33 
V. F. Krotov, Global Methods in Optimal Control Theory 
(Marcel Dekker, New York, 1996). [34 
A. Carlini, A. Hosoya, T. Koike, and Y. Okudaira, Phys. 
Rev. Lett. 96, 060503 (2006). [35 
A. T. Rezakhani, W.-J. Kuo, A. Hamma, D. A. Lidar, 
and P. Zanardi, Phys. Rev. Lett. 103, 080502 (2009). [36 
F. Verstraete, M. M. Wolf, and J. I. Cirac, Nat. Phys. 5, 
633 (2009). [37 
A. Sporl, T. Schulte-Herbruggen, S. J. Glaser, 
V. Bergholm, M. J. Storcz, J. Ferber, and F. K. Wil- [38 
helm, Phys. Rev. A 75, 012302 (2007). 

U. Schollwock, Rev. Mod. Phys. 77, 259 (2005). [39 
O. Romero-Isart and J. J. Garcia-Ripoll, Phys. Rev. A 
76, 052304 (2007). [40 
L. Jiang, A. M. Rey, O. Romero-Isart, J. J. Garcia-Ripoll, 
A. Sanpera, and M. D. Lukin, Phys. Rev. A 79, 022309 [41 
(2009). [42 
F. Galve, D. Zueco, G. Reuther, and S. Kohler, Eur. 
Phys. Journ. Special Topics 180, 237 (2010). [43 
C. DiFranco, M. Paternostro, and M. S. Kim, Phys. Rev. 
A 81, 022319 (2010). [44 
J.-S. Li, J. Ruths, and D. Stefanatos, J. Chem. Phys. 
131, 164110 (2009). [45 
J. Ruths and J.-S. Li (2011), arXiv:1102.3713. 

H. Rabitz, M. Hsieh, and C. Rosenthal, Science 303, 1998 [46 
(2004). 



H. Rabitz, T. S. Ho, M. Hsieh, R. Kosut, and M. Demi- 
ralp, Phys. Rev. A 74, 012721 (2006). 
T. Caneva, T. Calarco, R. Fazio, G. E. Santoro, and 
S. Montangero (2010), arXiv:1011.6634. 
T. Caneva, M. Murphy, T. Calarco, R. Fazio, S. Mon- 
tangero, V. Giovannetti, and G. E. Santoro, Phys. Rev. 
Lett. 103, 240501 (2009). 

W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. 
Flannery, Numerical recipes in C: the art of scientific 
computing (Cambridge University Press, 1992), 2nd ed. 
Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. 
Phys. 73, 357 (2001). 

G. Wendin and V. Shumeiko, in Handbook of Theoreti- 
cal and Computational Nanotechnology, edited by M. Ri- 
eth and W. Schommers (American Scientific Publishers, 
2006), p. Ch. 12, cond- mat/0508729. 

H. J. Lipkin, N. Meshkov, and A. J. Click, Nucl. Phys. 
62, 188 (1965). 

R. Botet and R. Jullien, Phys. Rev. B 28, 3955 (1983). 
A. Messiah, Quantum mechanics, vol. 2 (North-Holland, 
Amsterdam, 1962). 

S. Sachdev, Quantum Phase Transition (Cambridge Uni- 
versity Press, 1999). 

W. H. Zurek, U. Dorner, and P. Zoller, Phys. Rev. Lett. 
95, 105701 (2005). 

A. Polkovnikov and V. Gritsev, Nature Physics 4, 477 

(2008) . 

F. Pellegrini, S. Montangero, G. E. Santoro, and R. Fazio, 
Phys. Rev. B 77, 140404(R) (2008). 
S. Deng, G. Ortiz, and L. Viola, Europhys. Lett. 84, 
67008 (2008). 

C. D. Grandi and A. Polkovnikov (2009), 
arXiv:0910.2236. 

U. Divakaran, V. Mukherjee, A. Dutta, and D. Sen 

(2009) , arXiv:0908.4004. 

J. Dziarmaga, Adv. Phys. 59, 1063 (2010). 

A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalat- 

tore (2010), arXiv:1007.5331. 

T. Caneva, R. Fazio, and G. E. Santoro, Phys. Rev. B 
78, 104426 (2008). 

V. Giovannetti, S. Lloyd, and L. Maccone, Phys. Rev. A 
67, 052109 (2003). 

M. Murphy, S. Montangero, V. Giovannetti, and 

T. Calarco, Phys. Rev. A 82, 022318 (2010). 

V. Balachandran and J. Gong, Phys. Rev. A 77, 012303 

(2008). 



10 



[47] X. Wang, A. Bayat, S. G. Schirmer, and S. Bose, Phys. 

Rev. A 81, 032312 (2010). 
[48] F. Platzer, F. Mintert, and A. Buchleitner, Phys. Rev. 

Lett. 105, 020501 (2010). 
[49] J. I. Latorre, R. Orus, E. Rico, and J. Vidal, Phys. Rev. 

A 71, 064101 (2005). 
[50] P. Calabrese and J. Cardy, J. Stat. Mech. p. P04010 

(2005). 

[51] W. Dur, L. Hartmann, M. Hein, M. Lewenstein, and H.- 
J. Briegel, Phys. Rev. Lett. 94, 097203 (2005). 



[52] J. I. Latorre and A. Riera, J. Phys. A: Math. Theor. 42, 

504002 (2009). 
[53] C. Zener, Proc. Royal Soc. A 137, 696 (1932). 
[54] F. Caruso, S. Montangero, T. Calarco, S.F. Huelga, and 

M.B. Plenio. [arX"iv:1103.0929l 
[55] The generalization of the problem to the optimization 

of an overall unitary transformation is straightforward, 

averaging over the contributions of a complete set of basis 

of the Hilbert space H. 



